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ABSTRACT 


This  paper  analyzes  the  effect  of  optical  feedback  on  laser 
frequency  stability  and  the  acousto  optic  isolator  concept,  which 
was  demonstrated  previously.  The  frequency  puiling  effect  asso- 
ciated with  feedback  of  a signal  at  a frequency  other  than  the 
oscillator  frequency  has  been  shown  to  require  a nonlinearity 
such  as  saturation  in  the  laser  medium.  The  analysis  mathemati- 
cally corroborates  the  initial  acousto  optic  isolator  concept  and 
the  limited  experimental  data  available. 

In  the  study  of  the  acousto  optic  isolator,  it  was  deter- 
mined that  an  acceptable  analytic  expression  for  the  imaginary 
part  of  the  susceptibility  of  an  unsaturated  inhomogeneously 
broadened  line,  and  the  expression  for  the  susceptibility  of  a 
saturated  inhomogeneously  broadened  line,  were  not  available. 

The  expressions  which  were  desired  for  the  analysis  in  this  paper 
were  obtained  for  what  is  believed  to  be  the  first  time. 

The  results  of  this  paper  provide  a foundation  for  the 
acousto  optic  isolator.  Applicability  of  these  findings  shou-d 
increase  as  coherent  laser  radars,  which  require  isolation  to 
insure  accuracy  of  the  Doppler  information  in  the  return  signal, 
are  developed.  Additionally,  the  results  reported  herein  should 
prove  useful  irryview  of  the  increasing  need  for  stable  amplitude 
optically  pumped', submillimeter  sources. 
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SUMMARY 

This  paper  analyses  the  effect  of  optical  feedback  on  laser 
oscillator  frequency  stability  and  provides  a foundation  for  the 
acousto  optic  isolator.  This  should  find  increased  applicability 
as  coherent  laser  radars  are  developed  which  need  isolation  to 
insure  accuracy  of  the  Doppler  information  in  the  return  signal, 
and  also  as  requirements  develop  for  stable  amplitude  submilli- 
meter lasers  that  are  optically  pumped. 

Expressions  are  obtained  from  the  frequency  of  oscillation 
for  a laser  without  external  feedback,  and  with  feedback  of  a 
signal  at  a frequency  other  than  the  oscillator  frequency. 

These  expressions  are  directly  applicable  to  the  analysis  of 
the  acousto  optic  concept.  Both  homogeneously  broadened  lines 
and  inhomogeneously  broadened  lines  are  considered. 

The  frequency  of  oscillation  for  the  unsaturated  condition 
and  the  saturated  condition  is  determined.  Saturation  is  found 
to  change  the  oscillator  frequency  for  both  the  homogeneously 
and  inhomogeneously  broadened  lines. 

The  effect  of  feedback  on  the  oscillator  frequency  is  de- 
termined for  the  unsaturated  condition  and  found  to  not  in- 
fluence the  oscillator  frequency  because  of  the  lack  of  coupling 
between  the  feedback  signal  and  the  oscillator  frequency. 

For  the  homogeneously  broadened,  saturated  condition,  how- 
ever, the  frequency  of  oscillation  is  seen  to  be  pulled  by  the 
feedback  signal,  and  the  effect  is  proportional  to  the  intensity 
of  the  feedback  and  the  lineshape  of  the  transition  at  the  feed- 
back frequency.  This  condition  is  considered  to  be  a good  ap- 
proximation to  the  medium  in  a CO^  laser. 
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For  the  inhomogeneously  broadened  saturated  condition,  the 
expression  for  the  pulling  effect  can  be  obtained,  but  no  simple 
relation  between  the  frequency  of  oscillation  and  the  feedback 
characteristics  is  obvious. 

In  the  process  of  the  analysis  in  this  paper,  some  expres- 
sions were  derived  that  were  needed  to  determine  the  operating 
frequency.  These  include  the  complex  susceptibility  for  an 
unsaturated  inhomogeneously  broadened  line  and  the  complex  sus- 
ceptibility for  a saturated  inhomogeneously  broadened  line. 
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I.  INTRODUCTION 

After  the  usual  sources  of  laser  frequency  instability, 
such  as  power  supply  instabilities,  mechanical  vibrations, 
acoustical  effects,  temperature  variations,  and  discharge  in- 
stabilities have  been  reduced  to  a satisfactory  level,  there 
sometimes  remains  another  source  of  frequency  instabilit'  — 
feedback  of  laser  radiation  from  components  external  to  '■he 
laser  cavity.  This  feedback  can  degrade  the  performance  of 
systems  that  require  stable  frequency  laser  sources. 

In  some  laser  radar  applications,  a low  power  stable  laser 
oscillator  is  used  to  drive  a high  power  laser  amplifier  or  to 
injection  lock  a high  power-pulsed  laser  source.  The  low  power 
laser  provides  a frequency  stable  signal  which  is  especially  de- 
sirable in  systems  using  heterodyne  detection  to  extract  Doppler 
information,  or  simply  for  increased  sensitivity  in  detection. 

In  such  systems,  isolation  is  needed  between  the  stable  oscilla- 
tor and  the  high  power  device  to  prevent  unwanted  feedback  from 
the  high  power  source  to  the  stable  oscillator.  The  signal 
feedback  has  a tendency  to  pull  the  frequency  of  the  stable 
oscillator,  or  in  extreme  cases,  to  damage  the  oscillator. 

Optically  pumped  submillimeter  lasers  use  gases  that  have 
very  narrow  absorption  lines  at  the  typical  operating  pressure. 
Consequently,  the  pumps,  which  are  often  CO^  lasers,  must  be 
frequency  stable;  otherwise,  absorption  by  the  lasing  gas  varies 
and  the  submillimeter  signal  fluctuates  in  intensity.  Therefore, 
close  coupling  is  required  to  insure  a stable  submillimeter 
signal . 
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The  standard  technique  used  to  achieve  isolation  in  the 
microwave  or  optical  region  is  to  use  a nonreciprocal  device. 

The  design  of  most  isolators  in  these  regions  is  based  on  the 
principle  of  Faraday  rotation  (Refs.  1 and  2).  In  the  IR  region, 
it  is  usually  necessary  to  cool  a material  to  liquid  nitrogen 
or  liquid  helium  temperatures  in  order  to  observe  the  Faraday 
rotation  effect  (Ref.  3)«  As  a result,  the  Faraday  rotator  is 
not  used  extensively  in  the  IR  region.  The  currently  used  al- 
ternatives, however,  are  quite  unsatisfactory  (Ref.  4 ) , or  are 
of  limited  applicability  (Refs.  5 and  6).  These  include  using 
a polarizer  alone,  which  provides  isolation  from  stray  radia- 
tion of  the  opposite  polarization  only,  or  a quarter  wave  plate, 
which  generates  circularly  polarized  light  after  one  trip 
through  the  plate  and  linearly  polarized  light  of  the  opposite 
polarization  to  the  input  polarization  after  the  return  trip. 

The  acousto  optic  isolator  (AOI)  is  a concept  developed  to 
provide  nonreciprocal  isolation  between  stages  in  a master 
oscillator/power  amplifier  (MOPA)  chain.  The  concept  is  based 
on  the  fact  that  optical  or  IR  radiation  passing  through  an 
acousto  optic  device  is  partially  deflected,  and  the  frequency 
is  shifted  by  the  interaction  with  the  acoustical  wave.  If  the 
shifted  wave  is  returned  back  through  the  acousto  optic  device, 
the  wave  is  again  partially  deflected  and  the  frequency  shifted 
so  that  the  wave  has  a total  shift  of  twice  the  acoustical  fre- 
quency. This  effect,  which  has  been  demonstrated,  can  be  used 
to  prevent  frequency  pulling  of  a signal  by  spurious  radiation 
which  may  be  fed  back  in  a MOPA  chain. 

AOIs  may  become  important  in  future  laser  radar  systems  as 
the  output  power  increases  in  sophisticated  systems  which  must 
accurately  measure  the  frequency  content  of  the  return  signal, 
and  in  submillimeter  systems  requiring  stable  amplitude  trans- 
mitters which  are  optically  pumped. 
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Previous  work  in  this  area  has  qualitatively  predicted  the 
AOI  effect  and  experiments  have  been  performed  to  verify  the 
qualitative  predictions.  The  objectives  of  this  paper  are  to 
quantitatively  assess  the  effect  of  optical  feedback  on  laser 
frequency  stability  and  to  provide  a quantitative  analysis  of 
the  AOI.  The  quantitative  analysis  of  the  AOI  can  then  be  used 
to  predict  the  performance  of  a device  as  a function  of  input 
parameters . 

In  Section  II  the  acousto  optic  concept  is  explained  in 
more  detail.  In  Section  III-A,  the  general  equation  for  fre- 
quency of  an  oscillator  is  determined,  followed  by  a determina- 
tion of  the  actual  frequency  equations  for  homogeneously  and 
inhomogeneously  broadened  lines.  The  effects  of  saturation  are 
examined  in  Section  III-B. 

Finally,  the  effects  of  the  amplitude  and  the  frequency  of 
the  feedback  signal  on  the  oscillator  frequency  are  determined 
in  Section  III-C,  and  these  results  are  related  to  the  perform- 
ance of  an  AOI. 

Because  of  the  extensive  number  of  equations  involved,  the 
body  of  this  paper  essentially  summarizes  the  results  which  are 
derived  in  Appendix  A. 


II.  THE  ACOUSTO  OPTIC  ISOLATOR  CONCEPT 


The  AOI  concept  is  based  on  the  fact  that  an  acousto  optic 
modulator  shifts  the  frequency  of  the  radiation  that  is  deflected 
from  the  input  optical  beam  by  the  acousto  optic  interaction. 
Because  an  acousto  optic  modulator  has  180°  rotational  symmetry 
about  the  axis  corresponding  to  the  acoustic  beam  direction,  as 
illustrated  in  Fig.  1,  an  optical  rave  incident  from  the  output 
side  of  the  modulator  will  be  frequency  shifted  in  the  same  di- 
rection as  an  optical  wave  incident  on  the  input  side  of  the 
modulator.  Therefore,  radiation  reflected  back  through  the 
modulator  will  be  shifted  by  twice  the  modulation  frequency  of 
the  device,  as  indicated  in  Fig.  2. 

The  double  frequency  shift  can  be  easily  observed  by  homo- 
dyne detection  of  the  radiation  from  a gas  laser  that  has  been 
passed  twice  through  an  acousto  optic  modulator.  The  output  of 
the  homodyne  detection,  when  the  local  oscillator  radiation  and 
the  signal  to  be  detected  impinge  on  the  detection,  is  twice 
the  frequency  of  the  acoustical  input,  thus  verifying  the 
double  frequency  shift  (Ref.  7). 

Isolation  can  be  achieved  by  means  of  a modulator,  the 
acoustic  frequency  of  which  is  such  that  the  feedback  signal  is 
shifted  outside  of  the  gain  curve  of  the  oscillator  medium  as 
shown  in  Fig.  3-  In  the  case  of  a COp  oscillator,  which  has  a 
fluorescent  linewidth  Af  of  approximately  50  MHz,  the  reflected 
signal  can  be  shifted  outside  of  the  gain  curve  of  the  medium 
with  an  acoustical  signal  whose  frequency  is  approximately  30 
MHz.  In  this  frequency  range,  acoustical  signals  are  only 
slightly  attenuated  In  the  acousto  optic  medium  and  the  coupling 


- V" 


T"  PFSCEDIMS  PAGE  BLANK-NOT  FUMED 

c.  ••  «, 


7 


1 


I 


1 . 
i 


INPUT 

raoiation 


f 


m 


ACOUSTIC  WAVF 


■ 


f IS  LASER  FREQUENCY 
f IS  MODULATION  FREQUENCY! 

m i 


\m 

I 

ACOUSTIC  WAVE 


FIGURE  1.  Illustration  of  180°  Rotational  Symmetry  of  an 
Acousto  Optic  Material 
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FIGURE  2.  Illustration  of  Double  Frequency  Shift  of 
Optical  Beam  Due  to  Double  Pass  Through  an 
Acousto  Optic  Device 
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f - 2fm  f - fm  f 

11-2-76-3 

FIGURE  3.  Illustration  of  Shift  of  Feedback  Signal  of 

the  Gain  Curve  of  the  Oscillator.  For  a typical 
C02  osri nation  Af  = 50  MHz,  fm  = 30  MHz. 

An  extraneous  signal  that  enters  the  laser  medium  acts  as 
a loss  as  far  as  the  oscillator  signal  is  concerned  because  the 
extraneous  signal  depopulates  the  gain  medium.  This  loss  re- 
sults in  pulling  of  the  oscillator  frequency.  When  the  feed- 
back frequency  is  outside  of  the  gain  transition,  the  loss  due 
to  depopulation  is  reduced,  thereby  decreasing  the  frequency 
pulling,  or  increasing  the  stability  of  the  laser  oscillator. 
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III.  THEORY 

The  detailed  problem  to  be  analyzed  is  diagramed  in  Fig.  4 . 
When  a signal  other  than  one  at  the  frequency  of  oscillation  hits 
the  front  mirror  (assumed  to  be  the  only  port  of  entry  for  feed- 
back), some  of  the  radiation  is  reflected.  Also,  part  of  it 
enters  the  laser  cavity,  where  it  can  be  amplified  if  the  fre- 
quency of  the  radiation  is  within  the  gain  curve  of  the  laser 
medium.  The  amplification  occurs  at  the  expense  of  the  laser 
oscillating  frequency,  which  loses  available  energy  for  amplifica- 
tion stored  in  the  inverted  population.  This  effect  should  not 
be  as  pronounced  for  an  inhomogeneously  broadened  line,  since  there 
are  inverted  states  that  are  not  affected  by  the  laser  oscilla- 
tion frequency  in  the  inhomogeneously  broadened  case. 

The  problem  is  one  in  which  a cavity  exists  that  determines 
an  oscillation  frequency.  A driving  function  represented  by  the 
feedback  signal  is  input  to  the  cavity,  and  the  effect  of  the 
feedback  signal  whose  frequency  is  f „ , on  the  oscillator  frequency, 
f ^ , is  to  be  determined. 

To  solve  the  problem,  one  must  first  derive  the  general 
equation  for  the  oscillator  frequency  without  feedback,  and  then 
determine  the  frequency  for  the  type  of  transition  involved. 

The  oscillator  frequency  equation  with  feedback  is  derived  and 
compared  to  the  equation  without  feedback  to  determine  the  effect 
of  the  feedback.  This  equation  is  then  related  to  the  effect  that 
the  magnitude  and  frequency  of  the  input  AOI  signal  has  on  the 
optical  feedback  signal. 
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FIGURE  4. 


A.  UNPERTURBED  OSCILLATOR  FREQUENCY 

The  frequency  of  oscillation  of  a laser  has  been  approximately 
determined  by  considering  the  frequency  pulling,  caused  by  the 
phase  shift  per  transit  produced  by  the  refractive  laser  medium 
on  the  radiation  (Ref.  8),  or  by  considering  the  frequency  pulling 
due  to  dielectric  losses  which  are  equated  to  conductivity  cavity 
losses  for  a homogeneously  broadened  line  (Ref.  9). 

The  oscillating  frequency  can  be  determined  more  accurately 
by  solving  the  eigenvalue  equation.  It  is  derived  from  the  wave 
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equation  that  describes  the  electromagnetic  field  in  the  laser 
cavity,  by  substituting  the  material  equation  that  relates  the 
polarization  to  the  electric  field  into  the  wave  equation.  When 
this  procedure  is  followed,  as  shown  in  Appendix  A,  the  angular 
frequency  is  given  by 


2 . co 

w ) XIM  = t XRE 
c 


(2) 


+ 


where 

<oc  is  the  cavity  angular  frequency 
to  is  the  angular  frequency  of  oscillation 
t is  the  cavity  lifetime 

is  the  real  part  of  the  linear  susceptibility 

and 

is  the  imaginary  part  of  the  linear  susceptibility. 


This  equation  provides  the  general  solution  for  the  angular  fre- 
quency, to,  of  the  oscillator  whose  cavity  satisfies  Eq.  IB  in 
Appendix  A.  The  angular  frequency  depends  upon  both  the  real 
and  imaginary  parts  of  the  susceptibility.  (The  real  part  of 
the  susceptibility  of  a material  causes  a phase  shift  of  an  electro- 
magnetic wave  as  it  propagates  through  the  material.  The  fact 
that  the  phase  shift  is  frequency  dependent  results  in  the  dis- 
persion of  light  as  it  passes  through  the  material.  The  imaginary 
part  of  the  susceptibility  results  in  absorption  or  amplification 
of  waves.)  Therefore,  in  order  to  evaluate  this  equation  it  is 
necessary  to  determine  the  equation  for  the  susceptibility. 

Two  generic  types  of  transitions  occur  that  give  rise  to 
different  equations  for  susceptibility.  They  are  homogeneously 
broadened  and  inhomogeneously  broadened  transitions. 

^Appendix  A contains  the  entire  sequence  of  mathematical  com- 
putations for  this  paper,  from  which  specific  equations  have 
been  repeated  in  this  Section  for  the  convenience  of  the  reader. 
Therefore,  number  identification  of  equations  in  this  Section 
corresponds  to  the  placement  of  these  equations  in  Appendix  A. 
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1 . Homogeneous  Broadening 

When  the  macroscopic  polarization  of  a material  is  due  to 
the  sum  of  zhe  contributions  from  individual  atoms  which  have 
the  same  center  transition  frequency  and  linewidth,  the  absorp- 
tion or  spontaneous  emission  line  is  said  to  be  homogeneously 
broadened.  From  another  viewpoint,  the  effect  of  a group  of 
the  atoms  interacting  with  radiation  is  the  same  as  for  one  atom 
interacting  with  the  radiation  except  for  a multiplicative  con- 
stant. Homogeneous  broadening  results  whenever  the  relaxation 
processes  are  the  same  for  all  atoms.  These  relaxation  processes 
include  coupling  to  lattice  vibrations,  collisions,  and  spin- 
spin  coupling. 

The  values  of  XpiE  and  are  found  by  solving  the  dif- 

ferential equation  for  the  polarization  in  terms  of  the  electric 
field.  When  this  is  done  for  a homogeneously  broadened  line 
(Appendix  A),  the  results  are  as  shown  in  Fig.  5*  When  the  values 
for  Xpp  and  Xjm  are  substituted  into  Eq.  2,  the  general  equation 
for  the  oscillator  frequency,  then  the  equation  for  the  oscillat- 
ing frequency  is 

2 2 1/Tc  ,_2  2v 

“ ' “c  = 27TJ  w } (6) 

= T~A  u>(n2  - w2)  for  w «=  fi  , 

2 C 


which  shows  that  the  frequency  of  oscillation  assumes  a value 
between  the  center  of  the  transition  and  the  cavity  frequency. 

This  fact  is  illustrated  in  Fig.  6,  which  indicates  the  position 
of  the  oscillation  frequency  compared  to  the  spontaneous  emis- 
sion linewidth  and  the  cavity  linewidth. 

2 . Inhomogeneous  Broadening 

When  the  lineshape  of  a collection  of  atoms  or  molecules 
is  different  for  each  particle  because  of  slightly  different 
transition  frequencies  and/or  different  relaxation  times,  then  the 


I 


FIGURE  5.  Real  and  Imaginary  Parts  of  the  Complex 
Lorentzian  Lineshape  Factor 


line  is  said  to  be  inhoraogeneously  broadened.  Crystalline  in- 
homogeneities that  produce  different  local  fields  at  different 
lattice  points  in  a crystal,  nonuniform  magnetic  fields  envel- 
oping a spin  system,  or  the  Doppler  effect  due  to  motion  of 
atoms  or  molecules  in  a gas,  all  give  rise  to  inhomogeneous  broad- 
ening. The  resulting  lineshape  for  an  inhomogeneously  broadened 

'If  line  is  determined  by  summing  the  homogeneous  lineshapes  of  the 

individual  particles  involved  in  the  process.  The  result  of 
summing  a large  number  of  Lorentzian  lines,  as  would  be  expected 
by  the  central  limit  theorem,  is  a Gaussian  lineshape  (Ref.  10) 
(cf.  Fig.  7). 

The  same  procedure  is  followed  for  the  inhomogeneously 
broadened  case  as  was  used  for  the  homogeneously  broadened  case 
(Appendix  A)  to  obtain  the  susceptibility. 
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Absorption  line 


(a)  1/tc  « 2/T2  (b)  2/T2  « l/rc 

11-2-76-6 


FIGURE  6.  Frequency  Pulling  Effects  for  Cases  Where  Molecular  and 
Cavity  Bandwidths  Differ  Substantially.  The  case  shown 
in  (a)  is  in  general  applicable  to  the  laser,  whereas 
(b)  is  applicable  to  the  ammonia  cavity  laser  discussed 
in  Section  IV.  The  pulled  oscillation  frequency,  gj,  lies 
closest  to  the  narrower  of  the  two  lines.  Qj  and  Qc  are, 
respectively,  the  molecular  line  and  cavity  Q is  defined 
by  Qj  = ftT2/2,  Qc  e u^T.. 


gG’  8L 
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FIGURE  7.  Comparison  of  Gaussian  and  Lorentzian  Lines  of  the 
Same  Linewidth 
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A plot  of  the  real  and  imaginary  parts  of  the  susceptibility 
for  an  inhomogeneous ly  broadened  line  is  given  in  Fig.  8a  and  a 
comparison  of  the  real  susceptibility  for  a homogeneously  broad- 
ened line  and  an  inhomogeneously  broadened  line  is  made  in  Fig. 
8b. 


The  equation  for  the  frequency  of  oscillation  for  the  fre- 
quency of  an  inhomogeneously  broadened  line  is  obtained  when  these 
values  for  the  susceptibility  are  substituted  into  Eq.  2.  The 
equation  is  then 


(uj2-u )2)tt  = 2/tT 
c c 


[iUn2]1/2 


X 

e 


2 


dX 


(9) 


Although  this  equation  is  considerably  more  complicated  than 
Eq.  6,  it  nevertheless  shows  that  the  frequency  of  oscillation 
assumes  a value  between  the  center  of  the  transition  and  the 
cavity  frequency. 


B.  SATURATION  EFFECTS 

If  there  were  no  radiation  in  the  laser  cavity,  then  the 
population  inversion  would  be  determined  by  the  pumping  level 
and  the  relaxation  rates  for  transitions  from  the  upper  to  the 
lower  laser  level.  In  the  presence  of  electromagnetic  radiation, 
stimulated  transitions  occur  so  that  the  population  Inversion  in 
equilibrium  Is  reduced  from  the  value  obtained  without  radiation 
in  the  cavity.  As  the  level  of  radiation  increases,  the  value 
of  the  population  inversion  decreases.  For  a noninverted  popula- 
tion, energy  Is  absorbed  and  then  is  dissipated  in  relaxation 
processes-;  however,  for  a large  energy  absorption  the  system 
may  not  relax  to  equilibrium  as  fast  as  it  can  for  small  signals, 
so  the  system  saturates.  These  ideas  can  be  made  quantitative 
by  considering  the  differential  equation  for  the  population  dif- 
ference and  Including  the  effect  of  the  electric  field  on  this 
population  difference. 
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Complex  Susceptibility  of  an  an  Inhomogeneously  Broadened  Line. 

Inhomogeneously  Broadened  Line.  (See  Fig.  7 for  comparison  of  imaginary 

(See  Fig.  5 for  homogeneously  parts.) 

broadened  line.) 


I 


If  the  effects  of  the  electromagnetic  field  are  included, 
then  the  steady  state  solution  for  the  population  difference  is 
obtained  by  setting  the  time  dependent  parts  equal  to  zero.  The 
result  (Appendix  A)  is: 


N _N.  - ' 


T1 


co  e 

o 

~7sr  xtm 


(10) 


Equation  10  is  the  general  equation  for  the  steady  state  popula- 
tion difference.  It  shows  that  the  population  difference  is  a 
function  of  the  imaginary  part  of  the  susceptibility  of  the 
medium  and  the  field  intensity.  The  value  of  the  susceptibility, 
of  course,  is  itself  a function  of  the  population  difference. 

This  general  equation  can  be  used  to  determine  the  effects 
of  saturation  on  the  frequency  of  oscillation  by  substituting 
the  specific  expressions  for  the  imaginary  part  of  the  suscepti- 
bility, solving  for  the  population  difference,  and  substituting 
the  population  difference  into  Eq . 2 to  determine  the  oscil- 
lation frequency.  The  values  of  the  susceptibility  for  homo- 
geneously and  inhomogeneously  broadened  lines  can  now  be  sub- 
stituted . 

1 . Homogeneous  Broadening 

The  final  form  for  the  saturation  effect  depends  upon  the 
type  of  transition  involved,  since  the  form  of  the  suscepti- 
bility depends  upon  the  character  of  the  transition. 

The  effect  of  saturation  is  determined  by  substituting  the 
imaginary  part  of  the  susceptibility  for  the  specific  transition 
Into  Eq . 10  and  substituting  Eq . 10  Into  the  equation  for  the 
complex  susceptibility  of  the  transition.  The  saturation  effect 
obtained  by  this  procedure  for  a homogeneously  broadened  line 
(Appendix  A)  is  illustrated  in  Fig.  9. 
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11-2-76-11 

FIGURE  9.  Absorption  Constant  Under  Saturating  Conditions  for  a 
Homogeneously  Broadened  Line 


The  complex  susceptibility  is  affected  only  in  the  denomi 

p 

nator  by  the  change  in  |E|  or  to,  and  the  substitution  into  Eq 
2 gives 


Consequently,  the  equation  for  the  operating  frequency  is  un- 
changed explicitly  by  saturation  effects  on  the  susceptibility. 

A change,  however,  does  occur  in  the  cavity  frequency  oj^, 
since  the  cavity  frequency  is  related  to  the  index  of  refraction, 
which  depends  on  the  susceptibility.  When  this  factor  is  con- 
sidered, Eq.  2,  for  the  case  where  saturation  of  a homogeneously 
broadened  line  occurs,  becomes 


{[“co(1  - - “2j 
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^RE 


(15) 


Saturation,  therefore,  does  change  the  oscillator  frequency  from 
the  frequency  obtained  in  the  unsaturated  condition  for  a homo- 
geneously broadened  line. 

2 . Inhomogeneous  Broadening 

For  an  inhomogeneously  broadened  line,  the  equation  for  the 
population  difference  that  is  used  for  a homogeneously  broadened 
line  is  applicable  to  those  atoms  or  molecules  whose  transitions 
are  centered  at  a specified  frequency.  The  total  population 
difference  is  the  integral  of  the  specific  population  difference 
over  all  possible  frequencies.  Ordinarily,  when  the  integration 
is  performed  the  result  for  the  imaginary  part  of  the  suscepti- 
bility is  the  same  as  for  the  homogeneously  broadened  line,  ex- 
cept that  the  saturation  effect  depends  on  the  square  root  of 
intensity  rather  than  the  first  power.  In  order  to  calculate 
the  frequency  of  oscillation,  however,  the  real  term  is  needed 
as  well. 

If  the  same  procedure  as  for  the  homogeneously  broadened 
line  is  followed,  the  complex  susceptibility  Is  obtained  (Appendix 
A),  and  the  real  and  imaginary  parts  of  the  susceptibility  are 
plotted,  as  in  Figs.  10a  and  10b. 
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FIGURE  10a 


Real  Susceptibility.  Illustration  of  saturation  effect 
for  three  values  of  1/ 1 SAT  on  the  susceptibility  of  an 
inhomogeneously  broadened  line. 
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FIGURE  10b.  Imaginary  Susceptibility 

This  equation  for  the  susceptibility  is  considerably  more 
complicated  than  the  susceptibility  for  the  homogeneously  broad- 
ened case.  In  this  case,  the  effect  of  the  angular  frequency,  uj, 
and  the  intensity  effect  given  by  |E|  , are  different  for  the 
real  and  imaginary  parts  of  the  susceptibility.  Consequently 
the  frequency  equation 
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is  explicitly  affected 

by  the  saturation. 

(19) 


The  cavity  frequency,  w , is  also  affected  by  the  saturation 
for  the  inhomogeneously  broadened  line  as  given  by  Eq.  15: 
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so  that  saturation  has  a twofold  effect  on  the  oscillator  fre- 
quency for  an  inhomogeneously  broadened  line. 


C.  EFFECT  OF  FEEDBACK 

The  feedback  of  radiation  into  a laser  oscillator  can  be 
considered  as  the  insertion  of  an  electric  field  of  frequency 
other  than  that  of  the  oscillator,  as  indicated  in  Fig.  4.  This 
field  is  amplified  in  the  laser  medium,  thus  decreasing  the  popu- 
lation inversion  and  acting  as  a loss  in  the  medium.  It  is  im- 
portant to  recognize  that  the  loss  is  a dielectric  loss,  and  not 
a conductive  loss,  in  terms  of  the  macroscopic  wave  equation  that 
describes  the  processes.  This  is  due  to  the  fact  that  the  change 
in  population  causes  a change  in  x>  which  is  related  to  the  di- 
electric constant.  Since  the  amplification  of  the  feedback  depends 
upon  the  gain  profile,  off  resonance  radiation  should  not  be  as  ef- 
fective as  near  resonance  radiation  in  pulling  the  frequency. 

The  effect  of  feedback  on  the  oscillating  frequency  can  be 
analyzed  by  the  method  that  was  used  to  determine  the  oscillating 
frequency  of  a laser  without  external  feedback,  excot  that  now 
the  feedback  signal  must  be  included.  That  is,  the  differential 
equation  for  the  polarization  is  used  to  determine  the  suscepti- 
bility, and  this  result  is  substituted  into  the  cavity  equation 
to  obtain  the  canonical  equations  that  specify  the  frequency  of 
oscillation . 
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Unsaturated  Transitions 


1 . 


For  the  unsaturated  case,  the  analysis  shows  the  oscilla- 
tor frequency  equation  can  be  written  as 
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XIM  ^n^  = t XRE 
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(20) 


but  and  XRR  were  seen  to  be  functions  of  only  one  frequency, 

so  io^  is  Independent  of  ion . That  is,  when  the  medium  is  linear 
no  coupling  exists  between  the  two  frequencies;  therefore,  they 
act  independently.  Mathematically,  this  is  to  be  expected  since 
a linear  differential  equation  with  a sum  of  linear  inputs  has 
a solution  which  is  the  sum  of  the  solutions  due  to  the  individ- 
ual inputs;  however,  physically  it  is  somewhat  surprising  since 
the  inverted  population  is  changed  by  the  feedback  signal.  The 
dilemma  is  resolved  by  recognizing  that  the  change  in  population 
of  the  laser  medium  by  the  feedback  signal  requires  saturation. 
Otherwise,  the  population  supplying  energy  to  the  oscillator  fre- 
quency is  not  changed  by  the  feedback  signal. 


2.  Saturated  Transitions 


In  order  for  the  feedback  wave  to  interact  with  the  oscil- 
lator frequency,  in  general  a nonlinearity  is  required  in  the 
system.  This  nonlinearity  is  provided  by  the  saturation  effect. 
The  effect  of  the  feedback  signal  on  the  laser  oscillation  fre- 
quency can  be  determined  by  considering  the  equation  for  the 
population  difference,  which,  it  will  be  recalled,  is  applicable 
to  a homogeneous  line  or  to  a set  of  particles  in  an  inhomogeneous 
line  with  identical  energy  levels. 

Again,  the  analysis  shows  (Appendix  A)  that  the  change  in 
the  population  does  not  affect  the  frequency  of  oscillation 
through  the  explicit  susceptibility.  The  cavity  frequency, 
however,  does  change  so  that 
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where  is  the  frequency  of  oscillation  without  feedback.  From 

Eq.  25  the  oscillator  frequency  can  be  seen  to  be  linearly  re- 
lated to  the  value  of  the  transition  curve  at  the  feedback  fre- 
quency and  the  intensity  of  the  feedback.  Both  of  these  results 
have  been  experimentally  verified  as  indicated  by  Fig.  11. 
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Figure  11.  Comparison  of  Gain  Line  Profile 

and  Laser  Isolation  Achieved  versus 
Feedback  Offset  Frequency 
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For  the  inhomogeneously  broadened  line,  the  complex  suscep- 
tibility is  considerably  more  complicated,  as  can  be  seen  by  the 
analysis  in  Appendix  A. 

Fortunately,  for  the  case  of  a CO^  laser,  the  medium  is  not 
inhomogeneously  broadened.  This  is  due  to  the  fact  that  the 
Doppler  line  of  the  C02  laser  appears  to  be  homogeneously  broadened 
even  at  low  pressures  because  of  the  short  cross  relaxation  time 
of  the  molecules  across  the  velocity  profile  (Ref.  11).  The 
analysis  for  the  homogeneously  broadened  case  is,  therefore, 
more  likely  to  be  applicable. 
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APPENDIX  A 


DERIVATION  OF  EQUATIONS 

This  appendix  contains  the  mathematical  derivations  of  the 
equations  that  are  contained  in  the  body  of  this  paper.  The 
sections  in  this  appendix  correspond  to  the  sections  of  the  main 
text.  For  the  convenience  of  the  reader,  some  of  the  material 
from  the  main  text  is  repeated  below. 

III.  THEORY 

The  approach  to  solve  the  problem  diagramed  in  Fig.  4 
(p.  12)  requires  deriving  the  general  equation  for  the  oscil- 
lator frequency  without  feedback,  and  then  determining  the 
frequency  for  the  type  of  transition  involved.  The  oscillator 
frequency  equation  with  feedback  is  derived  and  compared  to  the 
equation  without  feedback  to  determine  the  effect  of  the  feed- 
back. 


A.  UNPERTURBED  OSCILLATOR  FREQUENCY 

The  frequency  of  oscillation  of  a laser  has  been  approxi- 
mately determined  by  considering  the  frequency  pulling,  caused 
by  the  phase  shift  per  transit  produced  by  the  refractive  laser 
medium  on  the  radiation,  or  by  considering  the  frequency  pulling 
due  to  dielectric  losses  which  are  equated  to  conductivity 
cavity  losses  for  a homogeneously  broadened  line.  (See  Refs. 

8 and  9 • ) 


acc 


The  laser  oscillating  frequency  can  be  determined  mor 
urately  than  in  Refs.  8 and  9 by  solving  the  eigenvalue 
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equation  which  is  derived  from  the  wave  equation  that  describes 
the  electromagnetic  field  in  the  laser  cavity.  The  equations  of 
interest  are  the  differential  equations  for  the  polarization, 
the  electric  field,  and  the  population  inversion.  When  these 
equations  are  derived,  as  in  Ref.  A-l,  the  results  are 


P + ft  P 


(N1  - N2)  E1 


E + — E+w2  E = - - PS 
tc  c e 


at  (Ni  + 


(N1  - N2)  - (Nx  - N2) 


— P 


eloc  ,(1C) 


where 


P is  the  polarization 

is  the  transverse  relaxation  time 

ft  is  the  angular  transition  frequency  of  the  medium 

h - h/2ir,  where  h is  Planck's  constant 

the  magnitude  of  the  induced  dipole  moment 

is  the  number  of  molecules/volume  in  the  lower  energy 
state 


is  the  number  of  molecules/volume  in  the  upper  energy 
state 


t or 

E is  the  local  electric  field 

E is  the  macroscopic  electric  field 

t is  the  cavity  lifetime 

u)c  is  the  cavity  angular  frequency 

e is  the  permittivity  of  the  medium 

Ps  is  the  polarization  for  the  transition  of  interest 

N^e  is  the  number  of  molecules/VOL  in  the  lower  energy 
state  in  equilibrium 

is  the  number  of  molecules/VOL  in  the  upper  energy 
state  of  equilibrium. 

If  the  polarization  is  expressed  in  terms  of  a linear 
susceptibility  by  the  equation 

P = X (<*>)  E 

(XrE  + 1 XIM)  D E 

where 

~ denotes  a complex  quantity 

is  the  permittivity  of  free  space 
Xr£  is  the  real  part  of  the  linear  susceptibility  x 
XIM  ls  the  imaginary  part  of  the  susceptibility, 

and  the  field  and  polarization  are  expressed  by 

— r (^t“icz) 

P = P e1  (<*>t-kz)  , 


^"Because  of  the  linearity  of  the  equations  involved,  the 
complex  form  for  the  electric  field  and  polarization  can 
be  used.  Later  in  the  paper  this  will  not  be  possibile. 


then  substituting  into  the  cavity  wave  Eq.  IB  for  the  electric 
field  gives 


(i^)2  + + “2  7 x(w)(id))2 

c c e 


+ i 

0) 

1 

a)  e 

o 

T=J 

X («*>) 

e 

+ i 

s-l 

(XRE 

1XiM) 

Since  the  right  side  of  this  equation  is  real,  the  imaginary 
part  of  the  left  side  must  be  zero,  i.e., 

(“c  “ “ ) XIM  = XRE  • (2) 

This  equation  provides  the  general  solution  for  the  angular  fre- 
quency, u>,  of  the  oscillator  whose  cavity  satisfies  Eq.  IB.  The 
angular  frequency  depends  upon  both  the  real  and  imaginary 
parts  of  the  susceptibility.  (The  real  part  of  the  suscepti- 
bility of  a material  causes  a phase  shift  of  an  electromagnetic 
wave  as  it  propagates  through  the  material.  The  fact  that  the 
phase  shift  is  frequency  dependent  results  in  the  dispersion  of 
light  as  it  passes  through  the  material.  The  imaginary  part  of 
the  susceptibility  results  in  absorption  or  amplification  of 
waves.)  In  order  to  evaluate  this  equation,  therefore,  it  is 
necessary  to  determine  the  equation  for  the  susceptibility. 

Two  generic  types  of  transitions  occur  that  give  rise  to 
different  equations  for  susceptibility.  They  are  homogeneously 
broadened  and  inhomogeneously  broadened  transitions. 


f 


1.  Homogeneous  Broadenin< 


When  the  macroscopic  polarization  of  a material  is  due  to 

the  sum  of  the  contributions  from  individual  atoms  which  have 
the  same  center  transition  frequency  and  linewldth-,  the  absorp- 
tion or  spontaneous  emission  line  is  said  to  be  homogeneously 
broadened.  From  another  viewpoint  the  effect  of  a group  of  the 
atoms  interacting  with  radiation  is  the  same  as  for  one  atom 
interacting  with  the  radiation  except  for  a multiplicative  con- 
stant. Homogeneous  broadening  results  whenever  the  relaxation 
processes  are  the  same  for  all  atoms.  These  relaxation  processes 
include  coupling  to  lattice  vibrations,  collisions,  and  spin- 
spin  coupling. 


The  polarization  for  the  homogeneously  broadened  line  is 
given  by 

" 0 " 0 on  I T p I 

P + s-P  + fi  P = =rL  — ( N - N-)  E , 

l 0 ft  D 1 c. 


(Nx  - IJ2)  E 


2 2 

where  L = [n  + 2/3]  is  the  Lorentz  correction  factor  for  an 
isotroptic  medium.  When  a solution  of  the  form 


E = E e 


i(iot  - kz) 


P = P e 


i(tot  - kz) 


is  assumed,  the  solution  to  the  differential  equation  is 

, ,2 


& L— — (N.  - N2) 
h 3 12 


(fi2  - uj2)  + i (2w/T2) 


where 


eQ  X U)  E , 


X(u>)  = L (Nx  - n2)  gL  (to,  a) 


and  g (uj,  Q)  is  the  complex  frequency  dependent  line  shape. 

L 

When  g (to,  ft)  is  written  in  terms  of  real  and  imaginary  parts, 

1j 


/ _ 2ft 

gL  (u,  ft)  = — 


(ft2-io2) 


i2ft 


2ui/T, 


11  ..(n2-u2)2  + (2w/T2)2  71  (ft2-w2)2  + (2u)/T2)? 


(4) 


so  that  the  linear  susceptibility  is 


X XRE  + 1 XIM 


where 


_ 1 T ^12  „ v 2ft(ft2-(o2) 

XRE-,TL—  <"l  V -(7-J5)2+(2u/T27 


-1 

LIM 


12 


(N, -N„ ) 


2ft2co/T, 


1 p ? ? ? 

1 2 (ft  -to  r + (2u>/T2r 


Ordinarily  the  assumption  of  near  reasonance,  to  «=  ft , is  made,  in 
which  case , 

~ 1 I u12 1 ^ 1 ~ 

p - 7i L ~^r~  (Ni  - V t*-»)  + '~itt2  e 

so 


~ , n,  i a - u 

gL  (u,  ft)  = - 


i 1 


l/T, 


11  (ft-u>)2  + (1/T2)2  n ( fi-ui ) 2 + (1/T2)2 


. (5) 


These  susceptibility  terms  corresponding  to  gL  (w,  ft)  are  plotted 
in  Fig.'  5 (p.  15).  When  the  values  for  xRE  and  XIM  are  sub_ 
stituted  into  Eq.  2,  the  general  equation  for  the  oscillator 
frequency,  then  the  equation  for  the  oscillating  frequency  is 

2 2 _ 1/tc  / 02  2, 

(o  - ioc  - 2/f~  - u)  ) 


= T2/t  c co(ft2  - oo2)  for  (o  ft>  ft  , 
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(6) 


which  shows  that  the  frequency  of  oscillation  assumes  a value 
between  the  center  oi  the  transition  and  the  cavity  frequency. 
This  fact  is  illustrated  in  Fig.  6 (p.  16),  which  indicates  the 
position  of  the  oscillation  frequency  compared  to  the  spontane- 
ous emission  linewidth  and  the  cavity  linewidth. 


2 . Inhomogeneous  Broadening 

When  the  lineshape  of  a collection  of  atoms  or  molecules 
is  different  for  each  particle  because  of  slightly  different 
transition  frequencies  and/or  different  relaxation  times,  then 
-the  line  is  saic  to  be  inhomogeneous ly  broadened.  Crystalline 
inhomogen6irt4j&s^that  produce  different  local  fields  at  dif- 
ferent lattice  point"s''±r*~-acrystal , nonuniform  magnetic  fields 
enveloping  a spin  system,  or  th£~'9qjojoler  effect  due  to  motion 
of  atoms  or  molecules  in  a gas  all  gi ve^i^i-ug^  to  inhomogeneous 
broadening.  The  resulting  lineshape  for  an  inhcThrs^eneously 
broadened  line  is  determined  by  summing  the  homogene ouS^ine- 
shapes  of  the  individual  particles  Involved  in  the  process* 

The  result  of  summing  a large  number  of  Lorentzian  lines,  as 
would  be  expected  by  the  central  limit  theorem,  is  a Gaussian 
lineshape  (Ref.  10). 


For  the  Doppler  broadened  case  where  the  number  of  molecules 
per  unit  volume  dM  in  a transition  frequency  range  of  dfl^  is  de- 
termined by  the  Maxwellian  thermal  equilibrium  velocity  distri- 
bution, dM  is 

dM  = Ny  <5q  (J^,  u)o ) df2±  , 


where  Ny  is  the  total  number  of  molecules  per  unit  volume  and 
g„  is  the  Gaussian  lineshape  factor.  The  equation  for  g_  is 

vJ  Cr 


gG  (ni’  wo^ 


where 


a)  is  the  transition  frequency  of  the  stationary  molecule  and 
o 

is  the  1/2  width  of  the  Gaussian  line. 

U 

The  Gaussian  line  is  compared  to  the  Lorentzian  line  in 
Fig.  7 (p.  16).  The  Lorentzian  line  is  seen  to  have  broader  skirts 
than  the  Gaussian  line,  but  a lower  peak  value,  so  that  the  energy 
is  more  concentrated  for  the  Gaussian  line. 

The  susceptibility  of  the  Doppler  broadened  line  is 
dP  = E dX  (w,  ^ ) , 

where  now  different  molecules  have  different  properties  accord- 
ing to  the  Maxwellian  distribution.  The  linear  susceptibility 
is,  therefore, 

| li,  p I ‘ ~ 

dX(u>,  = -^r  L — | — ^pll_p22^  gL^’  ^fii  ,U)o^dfii  ’ 


where 

Now  Ny  ( p 

dx  (w 
The  total 


is  the  probability  of  occupation  of  the  ith  state. 
, -,  -p' n 0 ) = N.-N_,  so  that 

11  dcL  1 d. 

, ,2 


, a, ) - 


i' 


susceptibility 


iWipl 

— ^ — (n1-n2)  gL(u,ni)gG(n1>  u>o)dni  . 

is 


X ( w ) = 


12 


(n1-n2) 


J 


Sp  ^ ^ 3 ^i  ) Sq  ^ > 


u>  )dfi, 
o i 


For  the  case  where  the  Lorentzian  lines  that  make  up  the  Gaussian 
line  are  much  narrower  than  the  overall  width  of  the  Gaussian 
lines,  gL(w,fr^)  can  be  approximated  by 
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I 


gL(u>,ni)  = ^ - i - 0>)  , 


so  that 


XIM(w)  ^ eo  L 3 ( N1  N2  ^ gG  ^w»  ' 


Prom  the  Kramer-Kronig  relations 

xRE(u,)  if  PP/  dQi  » 

— OO 

where  PP  is  the  principal  part  of  the  contour  integral.  The 
absorption  is  seen  to  have  the  same  value  as  for  the  homogeneously 
broadened  case,  except  for  the  lineshape  difference. 

The  real  part  of  the  susceptibility  is 

XIM 


Xre(“> 


- - ? pp/ 


n dfl, 

i 


r zjl 

- i PPJ  — 

~ _ m 


12 


(Nn  - 


N2)  gQ  (G1,u)q) 


-dfi , 


when  substitution  for  Xjm^)  ls  made.  Substituting  for 


gQ  u>Q)  yields 


v = _k_  ' 

xRElu'  - hea  3 


‘W  ppJ  S* 


P/ 


[CM/tr)g.n  2] 


1/2  -4S.n2(fi1-w 

e dtJ1 

i^-uJ  1 
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I 
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lu  I2  il/2 

(Nl-N?)KiiZ^]  PP 
heo  3 12  Au>g 


A-ll 
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P / e (AtV' 
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The  solution  for  the  integral  (Ref.  A-2) 


I = PP 


/ 


-4 ln2(n,-uj  )2 

i o 


Au)„ 


dn 


can  be  found  by  the  substitutions 
4«,n2 


Auk 


2 = C’  fii  " “'o  ~ S’ 


a)  - a)  = b > 
o 


which  simplify  the  integral  to 


I = PP  I e 


/ 


-CS 


S-b 


dS 


Now  the  denominator  Is 


1 _ S + b 
S-S  ' S2  - b3  ’ 


so  the  integral  can  be  written  as 


I 

I 


dS  • 
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Since  the  integral  of  an  odd  function  over  limits  that  are  sym- 
metrical about  the  origin  is  zero,  the  second  term  vanishes  so 
that 


I(C,b)  = 


c 

/ 


-cs 

b 7 dS 

S -b 


which  can  be  written  as 


-b2C  / e C^S  _b2) 

= be  I P p — dS 

S -b 


t)  c 

If  I(C,b)  is  multiplied  by  e , and  the  product  is  differentiated 
with  respect  to  C, 


d 


b2C 

{ I ( ,b)eD  °}_  d_ 
dC  ~ dC 


= -b 


dS 


The  integral  can  be  written  as  the  error  function  whose  value 
can  be  found  in  a table  of  definite  integrals.  The  integral 
is  equal  to 
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00 


-CX  ,Y  /ir 
e dx  = 

2/C 


so  that 


d tnc^bjeb  c).  _beCb2  It 


dC 


When  the  integration  is  performed 

b2C 


I(C,b)  e - 1(0, b)  = -b  /tt 
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Cb 


/C 


dC 


But  1(0, b)  = 0,  so 


I(C,b)  = -b  /tt  e 


-b2C 
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Cb 


/C 


dC 


Let 


2 2 
Cb  = X 


bcdC  = 2XdX 


so  that 


I(  C,b)  = - 2/?  e 


-b2C 


l 


b/C 


e x dX 


This  equation  can  be  identified  as  Dawson's  Integral,  which  is 
tabulated.  (See  Appendix  C and  Ref.  A-3. ) Notice  that  the 
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integrand  is  an  even  function,  so  the  integral  is  an  odd  function, 
i.e. , I(-b)  = -1(b). 


The  real  part  of  the  susceptibility  is,  therefore, 
2 


XRE(w)  = -r~ -P—  ( N -N„  ) 


heQ  3 


L ' M12  1 


[ ( 1} / 7T ) £n2] 
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1/2 
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Aw, 


I (c ,b ) 
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e dX  gr(u>,w  ), 


while  the  imaginary  part  has  been  given  by 

2 


xIM(w)  = _7T 


L |y12 


(N  -N  ) g (u),u  ) , 


he  3 '”1  2 6G  ' ’ o' 


(7B) 


where 


/ v2  4 Z,n2 

1 /P  ^ 2 

, , [(4/Ti)«,n2]1/£f„  ° Aw 

gG(“l"o>  ’ 4^ 6 


G . 


(8) 


A plot  of  the  real  and  imaginary  parts  of  the  susceptibility 
for  an  inhomogeneously  broadened  line  is  given  in  Fig.  8a  (p.  18), 
and  a comparison  of  the  real  susceptibility  for  a homogeneously 
broadened  line  and  an  inhomogeneously  broadened  line  is  given  in 
Fig.  8b  (p.  18). 

The  equation  for  the  frequency  of  oscillation  for  the  fre- 
quency of  an  Inhomogeneously  broadened  line  is  obtained  when  these 
values  for  the  susceptibility  are  substituted  into  Eq.  2.  The 
equation  is  then 
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Although  this  equation  is  considerably  more  complicated  than 
Eq.  6,  it  nevertheless  shows  that  the  frequency  of  oscillation 
assumes  a value  between  the  center  of  the  transition  and  the 
cavity  frequency. 


B.  SATURATION  EFFECTS 


If  there  were  no  radiation  in  the  laser  cavity,  then  the 
population  inversion  would  be  determined  by  the  pumping  level 
and  the  relaxation  rates  for  transitions  from  the  upper  to  the 
lower  laser  level.  In  the  presence  of  electromagnetic  radiation, 
stimulated  transitions  occur  so  that  the  population  inversion  in 
equilibrium  is  reduced  from  the  value  obtained  without  radiation 
in  the  cavity.  As  the  level  of  radiation  increases,  the  value  of 
the  population  Inversion  decreases.  For  a noninverted  population, 
energy  Is  absorbed  and  then  is  dissipated  in  relaxation  processes; 
however,  for  a large  energy  absorption  the  system  may  not  relax 
to  equilibrium  as  fast  as  it  can  for  small  signals,  so  the 
system  saturates.  These  ideas  can  be  made  quantitative  by 
considering  the  differential  equation  for  the  population  dif- 
ference 


ft'W  * 


(Nj-Nj)  - (Nj-Nj) 


The  first  term  is  the  power  delivered  to  the  dipole  medium,  the 
second  is  the  dissipation  of  the  power  to  the  surroundings , and 
the  right  side  of  the  equation  represents  power  delivered  by 
the  electromagnetic  field. 

When  no  field  Is  present,  the  steady  state  solution  is  ob- 
tained by  setting  the  first  term  equal  to  zero,  so  that 

(N1-N2)  = (N]_-N2)e  , 


4 


I 
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i.e.,  the  population  difference  is  the  equilibrium  value,  as 
is  expected.  In  thermal  equilibrium  this  is  the  Boltzmann 
distribution . 

If  the  effects  of  the  electromagnetic  field  are  included, 
then  the  steady  state  solution  for  the  population  difference  is 
obtained  by  setting  the  time  dependent  parts  equal  to  zero.  For 


n , Z i (wt-kz ) , + 

P = 1/2  P e + cc 


E = 1/2  E el(wt  kz)  + cc  , 


the  equation  becomes 
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(10) 


+At  this  point  it  is  necessary  to  use  the  complex  conjugate 
terms  because  of  the  terms  Involving  products  of  the  polar- 
ization and  the  electric  field.  Without  the  complex  con- 
jugates, no  steady  state  term  would  be  obtained. 
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Equation  10  is  the  general  equation  for  the  steady  state  popula- 
tion difference.  It  can  be  used  to  determine  the  effects  of 
saturation  on  the  frequency  of  oscillation  by  substituting  the 
specific  expressions  for  the  imaginary  part  of  the  suscepti- 
bility, solving  for  the  population  difference,  and  substituting 
the  population  difference  into  Eq.  2 to  determine  the  oscillation 
frequency.  The  values  of  the  susceptibility  for  homogeneously 
and  inhomogeneously  broadened  lines  can  now  be  substituted. 


1 . Homogeneous  Broadening 


The  final  form  for  the  saturation  effect  depends  upon  the 
type  of  transition  involved,  since  the  form  of  the  suscepti- 
bility depends  upon  the  character  of  the  transition.  The  imag- 
inary part  of  the  susceptibility  for  a homogeneously  broadened 
line  from  Eq . 3 is  found  to  be 


nL|n12r  (Nj— N2) 
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When  this  equation  is  substituted  into  Eq.  10,  the  result  is 


A-18 
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(11) 


This  equation  shows  that  the  change  In  population,  due  to  a 
change  in  gL,  occurs  in  the  denominator  of  the  equation.  Sub- 
stitution of  the  population  difference  from  Eq.  11  into  the 
complex  susceptibility,  which  is  given  in  Eq . 3 


X(w)  = 


ttL  | u12  | 2 (N1-N2 ) 


( tu  y ^2  ) 


yields 

X(w) 


irLly 


12 


(N1~N2} 


3£oA  SL  »LK,|2 

i + s Ti  eLlEl2 

3* 
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. (12) 


which  shows  that  both  terms  of  the  susceptibility,  the  imaginary 
part,  which  causes  the  absorption,  or  gain,  and  the  real  part, 
which  produces  a phase  shift,  are  dependent  on  the  frequency  of 
operation  co  and  the  strength  of  the  field  E.  The  saturation 
effect  for  a homogeneously  broadened  line  is  illustrated  in  Fig. 
9 (p.  20). 
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The  complex  susceptibility  is  affected  only  in  the  denomi- 

2 

nator  by  the  change  in  |E|  or  u>,  and  the  substitution  into  Eq . 

2 gives 


2. 
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Consequently , the  equation  for  the  operating  frequency  is  un- 
changed explicitly  by  saturation  effects  on  the  susceptibility. 

A change,  however,  does  occur  in  the  cavity  frequency  uc , 
which  is  given  by 


_ TTnc 

CO  = — = 

c nL 


(i*D 


where 


Since 


n is  the  longitudinal  mode  number 
c is  the  velocity  of  light 
n is  the  index  of  refraction  of  the  medium 
L is  the  resonator  length. 


n = l + x(w)  > 


the  cavity  frequency  is 
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for  small  values  of  Ax*  This  change  in  to  , due  to  the  field  de- 
pendent change  in  x>  results  in  a change  in  the  operating  fre- 
quency. Equation  2,  for  the  case  where  saturation  of  a homo- 
geneously broadened  line  occurs,  becomes 

,2  2,  _ (jo 

^ ^ XIM  tc  XRE 


'[to  (1 
I CO 


o 


(15) 


Saturation,  therefore,  does  change  the  oscillator  frequency  from 
the  frequency  obtained  in  the  unsaturated  condition  for  a homo- 
geneously broadened  line. 

2 . Inhomogeneous  Broadening 

For  an  inhomogeneously  broadened  line,  the  equation  for  the 
population  difference  that  is  used  for  a homogeneously  broadened 
line  is  applicable  to  those  atoms  or  molecules  whose  transitions 
are  centered  at  a specified  frequency.  The  total  population  dif- 
ference is  the  integral  of  the  specific  population  difference 
over  all  possible  frequencies.  Ordinarily,  when  the  integration 
is  performed  the  result  for  the  imaginary  part  of  the  suscepti- 
bility is  the  same  as  for  the  homogeneously  broadened  line  ex- 
cept that  the  saturation  effect  depends  on  the  square  root  of 
intensity  rather  than  the  first  power.  In  order  to  calculate 
the  frequency  of  oscillation,  however,  the  real  term  is  needed 
as  well. 

The  polarization  equation  for  the  set  of  particles  whose 
transitions  are  all  centered  at  the  same  frequency  is 
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where 
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for  a homogeneously  broadened  line.  The  total  effect  due  to 
all  particles  is  the  summation  over  each  set  of  particles  with 
different  center  transition  frequencies.  For  a continuous  dis- 
tribution of  transition  frequencies,  the  result  is 
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where  p(k)  is  the  probability  that  a particle  belongs  to  the 
set  whose  center  frequency  corresponds  to  k.  When  the  change 
in  population  due  to  the  saturation  effect  for  a homogenous 
line 
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is  substituted,  then  the  susceptibility  is 
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where 
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These  are  the  expressions  that  are  usually  used,  and  when  they 
are  substituted  into  the  expression  for  the  susceptibility 
Eq.  16,  it  becomes 
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over  the  region  of  integration  and  u>  & Q 
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The  evaluation  of  the  imaginary  part  of  the  susceptibility  in 
this  case  presents  no  problem,  so  little  difficulty  is  ordinarily 
encountered  since  the  imaginary  part  of  the  susceptibility, 
corresponding  to  the  absorption  or  amplification,  is  usually 
all  that  is  desired.  However,  the  real  part  of  the  suscepti- 
bility is  zero  under  these  assumptions  since  the  integrand  is 
odd.  The^pproach  of  assuming  p(ft,  ) is  a constant,  therefore, 
cannot  be  used 'Tor...  this  problem;  a more  rigorous  approach  is 
necessary . 

The  transitions  that  occur  ai>--aqiy  frequency  can  be  con- 
sidered to  be  the  sum  of  the  transit ions'S^at  can  occur  at  that 
frequency  due  to  all  of  the  homogeneously  broaTiejqed  lines  that 
make  up  the  inhomogeneously  broadened  curve,  weighted-'-tiy  the 
shape  of  the  inhomogeneously  broadened  curve.  The  weighting' — 
function  for  this  case  is 
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The  equation  for  the  susceptibility  is,  then, 
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It  would  appear  that  both  integrals  could  be  evaluated  by 
a contour  integration  where  the  contour  includes  the  real  axis 
and  the  upper  half  plane.  The  exponential  term,  however, 
diverges  for  large  imaginary  values  so  that  a reasonable  con- 
tour is  difficult  to  find. 

Various  attempts  to  alter  the  equation  to  fit  the  tabulated 
integrals  failed.  A complicated  solution  using  a power  series 
expansion  was  found.  The  complexity  encouraged  search  for  an- 
other solution,  which  was  finally  successful. 

The  solution  is  to  make  a substitution  which  causes  the 
integrals  to  take  a form  given  in  Ref.  11. 

The  first  integral 
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so  that  the  first  integral  becomes 
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are  substituted,  the  integral  is 
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The  second  Integral  is  found  in  a similar  way 
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Let 
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The  equation  for  the  susceptibility  is,  therefore, 
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The  real  and  imaginary  parts  of  Eq.  18  are  plotted  in  Figs.  10a 
and  10b  (pps.  22  and  23). 

This  equation  for  the  susceptibility  is  considerably  more 
complicated  than  the  susceptibility  for  the  homogeneously  broad- 
ened case.  In  this  case  the  effect  of  the  angular  frequency, 

(*),  and  the  intensity  effect  given  by  | E ( , are  different  for  the 
real  and  imaginary  parts  of  the  susceptibility.  Consequently 
the  frequency  equation 
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is  explicitly  affected  by  the  saturation. 
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since 


The  cavity  frequency,  u)fi , is  also  affected  by  the  saturation 
for  the  inhomogeneously  broadened  line  as  given  by  Eq.  15 
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so  that  saturation  has  a twofold  effect  on  the  oscillator  fre- 
quency for  an  inhomogeneously  broadened  line. 


C.  EFFECT  OF  FEEDBACK 


The  feedback  of  radiation  into  a laser  oscillator  can  be 
considered  as  the  insertion  of  an  electric  field  of  frequency 
other  than  that  of  the  oscillator,  as  indicated  in  Fig.  4 (p. 

12).  The  effect  of  feedback  on  the  oscillating  frequency  can  be 
analyzed  by  the  method  that  was  used  to  determine  the  oscillating 
frequency  of  a laser  without  external  feedback,  except  that  now 
the  feedback  signal  must  be  included.  That  is,  the  differential 
equation  for  the  polarization  is  used  to  determine  the  suscepti- 
bility, and  this  result  is  substituted  into  the  cavity  equation 
to  obtain  the  canonical  equations  that  specify  the  frequency  of 
oscillation. 
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Coefficients  of  the  same  frequency  can  be  equated  because  of  the 
orthogonality  of  the  two  input  frequencies. 


The  results  are 
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When  the  polarization  and  electric  field  are  substituted  into 
the  differential  equation  for  the  field  in  the  cavity 
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the  result  is 
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Again,  the  terms  of  the  same  frequency  can  be  equated  and  the 
expression  for  the  susceptibility  can  be  substituted  so  that  the 
canonical  equations  become 
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Each  of  these  equations  is  identical  in  form  to  the  oscillator 
frequency  equation  so  that  the  frequency  equation  can  be  written 
as 
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but  Xjjv-j  and  were  seen  to  be  functions  of  only  one  frequency, 
so  is  independent  of  • That  is,  when  the  medium  is  linear 
no  coupling  exists  between  the  two  frequencies;  therefore  they 
act  independently.  Mathematically,  this  is  to  be  expected  since 
a linear  differential  equation  with  a sum  of  linear  inputs  has 
a solution  which  is  the  sum  of  the  solutions  due  to  the  individ- 
ual inputs;  however,  physically  it  is  somewhat  surprising  since 
the  inverted  population  is  changed  by  the  feedback  signal.  The 
dilemma  is  resolved  by  recognizing  that  the  change  in  population 
of  the  laser  medium  by  the  feedback  signal  requires  saturation. 
Otherwise,  the  population  supplying  energy  to  the  oscillator 
frequency  is  not  changed  by  the  feedback  signal. 

2 . Saturated  Transitions 

In  order  for  the  feedback  wave  to  interact  with  the  oscil- 
lator frequency,  in  general  a nonlinearity  is  required  in  the 
system.  This  nonlinearity  is  provided  by  the  saturation  effect. 
The  effect  of  the  feedback  signal  on  the  laser  oscillation  fre- 
quency can  be  determined  by  considering  the  equation  for  the 
population  difference 
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where  P and  E each  consist  of  two  frequencies.  This  equation, 
it  will  be  recalled,  is  applicable  to  a homogeneous  line  or  to 
a set  of  particles  in  an  inhomogeneous  line  with  identical  energy 
levels . 

The  polarization  and  electric  field  are 
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The  time  independent  terms  are 
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The  steady  state  or  time  independent  part  of  the  equation  for 
the  population  difference  is,  therefore. 


(N1-N2)-(N1-N2)' 


2 n 


f~*  ~ \ /~  ~ \ 

l*El-PrElj+  u)2^P2*E2-P2*E2/ 


Now 


A-33 


^ : — ' 


a 


P = ?!  + p2  = X1E1  + *2E; 


(N1-N?)-(N1-N2)< 

tT 


2 htt  V 


'X1E1E1/ 


iX2E2E2*x2E2E2 


f,(l  U)lxlM^a)l^Ell  +w2XIM('a,2^E2  I 


where  all  time  dependent  terras.  Including  the  low  frequency  com- 
ponents generated  by  the  mixing  of  two  optical  frequencies,  have 
been  ignored.  Substitution  of  the  value  of  xIM(w  ) from  Eqs . 3 
and  5 , 
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The  new  population  difference  is 
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where  I = C | E | is  the  intensity  of  the  field  at  frequency 
co^.  This  equation  can  be  seen  to  be  a generalization  of  Eq.  11 
for  the  case  where  two  signals  are  influencing  the  oscillator. 
If,  therefore,  saturation  has  occurred  without  feedback,  as  it 
does  in  gas  lasers,  then  the  population  is  determined  by 
gLl(sL2  0) ’ so  the  frequency  is  fixed.  Now  feedback  by  an- 
other signal  changes  the  population  due  to  g . If  the  satura- 
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tion  term  is  small,  then 
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for  u2  , 


If 

where  Kw  = T2/tt,  which  shows  that  the  value  of  N -N2  is 

linearly  related  to  I2>  the  intensity  of  the  feedback  signal, 
and  (io2,n),  the  lineshape  of  the  transition.  When  the  new 
population  difference  is  substituted  into  Eq.  3,  the  expres- 
sion for  the  complex  susceptibility  becomes 
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The  exact  expression  for  the  frequency  pulling  or  pushing 
of  the  oscillator  due  to  the  feedback  is  obtained  by  substituting 
the  change  in  susceptibility  due  to  the  population  change  into 
the  general  equation 
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Again,  the  change  in  the  population  does  not  affect  the 
frequency  of  oscillation  through  the  explicit  susceptibility 
expression  since,  from  the  general  frequency  equation, 
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and  the  ratio  XRE/XIM  cancels  any  change  in  the  real  and  imaginary 
parts  of  the  susceptibility  due  to  population  changes.  The  cavity 
frequency,  however,  does  change  so  that 
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which  for  small  changes  in  the  frequency  from  the  cavity  frequency 
is  given  by 
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which,  from  Eq.  22,  becomes 
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where  is  the  frequency  of  oscillation  without  feedback.  Prom 

Eq,  25  the  oscillator  frequency  can  be  seen  to  be  linearly  re- 
lated to  the  value  of  the  transition  curve  at  the  feedback  fre- 
quency and  the  intensity  of  the  feedback.  Both  of  these  results 
have  been  experimentally  verified  as  indicated  by  Pig.  11  (p. 

26). 

For  the  inhomogeneously  broadened  line,  the  complex  suscep- 
tibility is  considerably  more  complicated.  In  this  case  the 
complex  susceptibility  is  found  by  substituting  Eq.  21  into  Eq. 
16,  so  that 
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if  w ~ then  the  values  of  the  denominators  are  such  that 
(fii  - w1)2  + (1/T2)2  m (fii  - .2)2  + (1/T2)2 


for  u>^  ,uj0  ~ . 

The  equation  for  the  complex  susceptibility  can  then  be 
simplified  to 

2 


ti  Ll  U12  1 

x(“i>  = W~  T <W 


ek 


/ 


1 (nk  wl) 


p(fik)  dk 


(Pk-Wl)2  + (1/T2)2  + ^ ~ (u)1|E1|2  + u)2|E2|2) 

2 1c 


0 

/ 

J — o 


1/T, 


p(^)dk 


I^)2  + a/T2)2  + (Wi|Ei,2  + w2 |e2| 2) 
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If  fi.  in  the  right  side  of  the  denominator  can  be  assumed  to  be 
constant  over  the  range  of  integration,  then  x(w1)  is  the  same  form 
as  Eq.  17B,  where 

/ 1 \2  L|  u,  5 1 T-ji  P P 

* - (tt)  + — r-  r;  (“i i Ei'  + “2iE2n  > 

so  that  the  final  solution  is  Eq.  19,  with  y as  given  above. 

This  equation  is  difficult  to  evaluate  compared  to  the  equation 
for  the  homogeneously  broadened  case. 

Fortunately,  for  the  case  of  a C0o  laser,  the  medium  is  not 
inhomogeneously  broadened,  so  the  analysis  for  the  homogeneously 
broadened  case  is  more  likely  to  be  applicable. 
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APPENDIX  B 

UNSUCCESSFUL  ATTEMPTS  TO  FIND  I 
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UNSUCCESSFUL  ATTEMPTS  TO  FIND  I 


The  first  approach  toward  solving  this  equation  is  to  search 
for  a definite  integral  in  a table.  After  no  satisfactory  tab- 
ular form  is  found,  a straightforward  contour  integral  is  at- 
tempted. 


The  integral 


I = PP  I e 


/ 


oo  -A2n2  ( ) ‘ 


(Aojp)' 


is  equal  to 


I = lim  J e 

R-»oo 


■I 

-R 


D -4dn2 (ft . —to  )‘ 
R i o 


Aw, 


This  integral  has  a pole  on  the  positive  real  axis  because  w, 
the  angular  frequency,  is  positive  and  real.  Since  an  optical 
frequency  is  involved,  it  is  also  apparent  that  u>>>0.  The  func- 
tion can  be  visualized,  as  indicated  in  Fig.  B-l,  as  the  product 

of  a Gaussian  Line  centered  at  a fixed  value  w , and  the  term 

o 

1 /(ft.-ui)  where  w is  variable  but  close  to  w . Because  the  in- 
l o 

tegral  has  a real  pole,  the  theory  of  residues  cannot  be  used 
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since  the  denominator  of  an  integral  to  be  evaluated  by  the 
theory  of  residues  cannot  have  a real  zero  of  any  order 
(Ref.  B-l ) . 


FIGURE  B-l.  Illustration  of  Two  Functions  of  Which  Integrand  is 
Composed.  Because  Gaussian  function  is  symmetrically 
about  u*),  and  (1/fi-j-w)  is  anti symmetri cal  about  co0, 
the  integral  vanishes  when  oj  = oj  . 


Rather  than  a contour  integral  which  includes  the  pole,  a 
contour  integral  excluding  the  pole  can  be  attempted.  If  a 
contour  is  chosen  that  excludes  the  pole,  then  the  integral  Is 
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-4Sln2(Reie-wo)2 

7~ie 
Re  -a) 


Rd0  , 


which  becomes 

(Jl) 
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-4&n2  ( £7^— toQ ) ‘ 


1 ■ 1 e A“a 


o - + 

. -aj  i 


/ 


-4£n2(f2^-a)o)‘ 

7 2 

Aw„ 


-7: — — dfi, 

ft,  —10  i 


CO 


= 0 


By  inspection  of  Fig.  B-l,  it  is  clear  that  the  integral  is  zero 
for  <*)  = uj  , but  must  have  a value  for  co  ¥■  u)Q.  Also,  the  value 
of  the  integral  is  antisymmetrical  for  values  of  u about 

Another  approach  is  to  change  the  form  of  the  integral  so 
that  it  takes  the  form  of  a tabulated  integral.  If  the  substi- 
tution X = - a)  is  made 
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R 

-4^n2(ni-wo)2/Aa)G2 

I = llm  / £ dfi 

R-*-®  | n1-u> 


L 


R-k) 

f 

= lim  I e_ 
1 R-*-»  I 

-R-k) 


-4«.n2(X+k)-k)  r/Ak)p2 
o u 


dX 


R-k) 


= lim  # e 

R+oo 

’’’-R-k) 


-4«.n2[X2+2X(k)-k)Q)  + (k)-k)Q)2]/Ak)G2 


dX  * 


which  is  not  available  in  tables  either. 
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APPENDIX  C 

TABLE  OF  DAWSON'S  INTEGRAL 


This  Table  is  reproduced  from  Ref.  11  for  the  convenience 
of  the  reader. 
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ERROR  FUNCTION  AND  FRESNEL  INTEGRALS 
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0. 

59588 

1008 

2 

0.  08 
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2 
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0. 
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0. 
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0. 
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0. 

35178 

01580 

0.  055 

0. 
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4 
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0. 
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0. 
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0.  035 

0. 
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1.  90 

0. 
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0. 
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6 

0.92 
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1.92 

0. 

31847 

19293 

0.  920 

0. 
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7 

0.94 

0.  54090 
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1.  94 

0. 

31405 
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8 

0.  96 

0.  54036 
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1.  96 

0. 

30973 
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0. 
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50580 

1.  98 

0. 
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0. 
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00 
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APPENDIX  D 

TABLE  OF  ERROR  FUNCTION  FOR  COMPLEX  ARGUMENTS 

This  Table  Is  reproduced  from  Ref.  11  for  the  convenience 
of  the  reader. 


Tahir 


ERROR  FUNCTION  FOR  COMPLEX  ARGUMENTS 
«<*)  c !,erfc(  iz)  i *+iy 


*»c(z) 

<»4  (Z) 

*U'{Z) 

&U'(  z I 

' '»<’(*) 

fti'(z) 

<nU) 

&\v(  Z 1 

flC(Z) 

> 

» 

0 

X 

0.1 

X 

0.2 

X 

0.3 

X 

0.4 

0.0 

1.000000 

0.000000 

0.990050 

0.112089 

0.960789 

0.219753 

0.91)9)1 

0.318916 

0.852144 

0.40615) 

0.1 

0.896457 

0.000000 

0.888479 

0.094)32 

0.86498) 

0.185252 

0.827?4b 

0.269600 

0.777267 

0.344680 

0.? 

0.809020 

0.000000 

0.802667 

0.080079 

0.78)8)8 

0.15740) 

0.752895 

0.22965) 

0.712146 

0.29465) 

0.3 

0.734599 

0.000000 

0.729)37 

0.068410 

0.71)801 

0.1)47)9 

0.688720 

0.1970)7 

0.65S244 

0.25)bl  ) 

0.4 

0.670788 

0.000000 

0.66646) 

0.058897 

0.65)680 

0.116147 

0.6)2996 

0.17020) 

0.605295 

0.219706 

0.5 

0.615690 

0.000000 

0.612109 

0.051048 

0.60181) 

0.100782 

0.584))) 

0.147965 

0.56125? 

0.191500 

O.b 

0 567805 

0.000000 

0.564818 

0.044524 

0.555974 

0.08799) 

0.541605 

0.129408 

0.522246 

0.167880 

0.7 

0.525930 

0.000000 

0.57)42) 

0.0)9064 

0.518991 

0.07  7275 

0.50)896 

0.11)821 

0.407556 

0.14  7975 

0.8 

0.489101 

0.000000 

0.486982 

0.0)4465 

0.480697 

0.0682)5 

0.4  7045? 

0.10064  7 

0.456579 

0.1)1101 

0.9 

0.456532 

0.000000 

0.4547)1 

0.0)0566 

0.449)83 

0.0608b) 

0.440655 

0.089444 

0.428808 

0.116714 

1.0 

0.427584 

0.000000 

0.426044 

0.027242 

0.421468 

0.054014 

0.41  3989 

0.079864 

0.40)818 

0.104)80 

1.1 

0.401  730 

0.000000 

0.400406 

0.024)9? 

0.3964  70 

0.048)9) 

0.390028 

0.071628 

0.381250 

0.09)752 

1.? 

0.3785)7 

0.000000 

0.377)9) 

0.0219)4 

0 373989 

0.04354? 

0.3684]? 
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0.36C799 

0.084647 

1.3 

0.35764) 

0.000000 

0.356649 

0.019805 

0.353691 

0.0)9)36 

0.3488)9 
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0.342206 

0.076538 

1.4 

0.3)8744 

0.000000 

0.337876 

0.01  7951 

0.335294 

0.035671 

0.3)1054 

0.0529)6 

0.325248 

0.069538 

1.5 

0.321585 

0.000000 

0.320825 

0.016329 

0.318861 

0.0)246) 

0.3148)9 

0.048210 

0.309736 

0.063)9) 

1.6 

0.30595) 

0.000000 

0.  306284 

0.014905 

0.30)290 

0.02964) 

0.300009 

0.044051 

0.295506 

0.057978 

1.7 

0.79166) 

0.000000 

0.291072 

0.01)648 

0.289309 

0.02  7154 

0.286406 

0.040)77 

0.28241  7 

0.05)186 

1.8 

0.278560 

0.000000 

0.2780)5 

0.0128)6 

0.276470 

0.024948 

0.27)892 

0.0)7118 

0.270346 

0.048931 

1.9 

0.266509 

0.000000 

0.266042 

0.011847 

0.264648 

0.022987 

0.262350 

0.0)4?]7 

0.259186 

0.0451 39 

2.0 

0.255)96 

0.000000 

0 264978 

0.010664 

0.253732 

0.0212)6 

0.251677 

0.031626 

0.248844 

0.041  748 

2.1 

0.745119 

0.000000 

0.244  746 

0.009874 

0.24)628 

0.019669 

0.241  78) 

0.029)04 

0.2)92)9 

0.0)8706 

2.2 

0.2)559) 

0.000000 

0.23528b 

0.009165 

0.2)4251 

0.018260 

0.23259? 

0.027217 

0.230)00 

0.0)5968 

2.3 

0.226742 

0.000000 

0.226438 

0.008826 

0.2258)1 

0.016991 

0.22405) 

0.025)35 

0.22196) 

0.0)3498 

2.4 

0.218499 

0.000000 

0.218724 

0.007949 
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0.015848 
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0.02)63) 

0.214172 

0.0)126) 

2.5 

0.210806 

0.000000 

0.210657 

0.007427 

0.209813 

0.014006 

0.208582 

0.022090 

0.206879 

0.0292  34 

2.6 

0.20)61) 

0.000000 

0.20)387 

0.00696? 

0.202710 

0.01)86? 

0.201589 

0.020687 

0.2000)9 

0.027)89 

2.7 

0.196874 

0.000000 
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0.006520 
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0.01)002 

0.195028 

0.019409 

0.19)61) 

0.025706 

2.8 

0.190^49 

0.000000 

0.190360 

0.006126 

0.189796 

0.012216 

0.188861 

0.018241 

0.187566 

0.024168 

2.9 

0.18460? 

0.000000 

0.184429 

0.005764 

0.18)912 

0.011498 

0.18)056 

0.01717? 

0.181868 

0.022759 

3.0 

0.179001 

0.000000 

0.178842 

0.0054)'' 

0.178)68 

0.0108)9 

0.177581 

0.016192 

0.1  76491 

0.021466 

\ 

X 

0.5 

X 

0.6 

X 

0.7 

X 

0.8 

1 

0.9 

0.0 

0.778801 

0.478926 

0.697676 
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( 0.5124242 

0.05176536  \ 

10  " 

\z?  0.2752551* 

*•’- 2.724745/' 

no) 

D-6 


I 


< mil'l  l A /.KUOS  Ol 

mi  i mu  in  1 1 muon 

1 ilili  7. Ill 
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2.  83974  105  3.  17562  810 

8 

4.  84797  031 

5.  10158  804 

4 

3.  33546  074  3.  64617  438 

9 
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5.  68883  744 
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erf  ( ',.)*(> 

"•  1 1,  >\ 

ln(’\ 2(' 
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-o.;2A  V"~  2}* 
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2/ 

From 

H K.  Salzer,  Complex  zeros  of  the  error  function,  .1  Franklin 

Inst.  200,  209  211 

1955 

(with  permission). 
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MAXIMA  AM)  MIMMA  OF  FRFSNFI*  IINTKOR ALS  Tabl«-7.I2 


Mn  C(V4»  ' 1) 

**1„  r(V4"  • 3) 

■>/;  ■ N( V4"  + 2)  *>*;  .S :(yf4 

1"  • 4 ) 

n 

l/„ 

"'n 

Mn 

r**f, 

0 

0. 779893 

0.  321056 

0.  713972 

0.  343415 

1 

1 1 

0.  640807 

0.  380389 

0.  628940 

0.  387969 

L*- 

2 

0.  605721 

0.  404260 

0.  600361 

0.  408301 

3 

0.  588128 

0.  417922 

0.  584942 

0.  420516 

4 

0.  577121 

0.  427036 

0.574957 

0. 428877 

5 

0.  569413 

0.  433666 

0.  567822 

0.  435059 

■ 

2 

*7(4*1  *1)J  3 
*3(4*1  3 I)'* 

1 *2(4*h3)2*3 
m"2  *3(4n3‘3)'f2 

1 »J(4n  * 2)2  3 . 1 *1)’  3 

"”'2  + »'(4'i  m"~2  32* ’(n  * 1)V* 


From  G N.  Watson,  A treatise  on  the  theory  of  Bessel  functions,  2d  ed.  Cambridge  Univ. 
Press,  Cambridge,  England,  1958  (with  permission). 
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